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We investigate the scalar Green function for spherically symmetric spacetimes expressed as a 
coordinate series expansion in the separation of the points. We calculate the series expansion of 
the function V{x,x') appearing in the Hadamard parametrix of the scalar Green function to very 
high order. This expansion is then used to investigate the convergence properties of the series and 
to estimate its radius of convergence. Using the method of Pade approximants, we show that the 
series can be extended beyond its radius of convergence to within a short distance of the normal 
neighborhood boundary. 



I. INTRODUCTION 



Quasi-local series expansions - expansions in the separation of two points x and x' - are a frequently used tool for 
calculations of fields on curved spacetimes. Often, as a final step in the calculation, the coincidence limit, x' — *■ x, 
is taken. In these cases, the precise convergence properties of the series is of little interest. However, there are cases 
where we would like the points to remain separated [TJ HJ |3J 31 [5] . In particular, we are motivated by the calculation 
of the quasilocal (QL) contribution to the scalar self-force [6l|7] (for a review, see [8l|9]) on a scalar particle. 



/QL(^(T)) = hmg2 / V-Gret{z{T),z{T'))dT', (1.1) 

'r-Ar 



where z(t) describes the worldline of the particle. This requires a quasi-local expansion of the retarded Green function, 
Gret{x, x'), which is a solution of the scalar wave equation with point source, 

(□-gi?)G',et(x,xO = -4^ ^ V- ' (1-2) 

where ^ is the curvature coupling constant and R is the Ricci scalar. As expression ( |1.1| requires the Green function 
for the points separated up to an amount At along a world-line, it begs the question: how large can the separation 
of the points be before the series expansion is no longer a valid representation of the Green function? 

To the authors' knowledge, this question has not yet been quantitatively answered. It is well known that the 
Hadamard parametrix for the Green function (upon which quasi-local calculations are based) is valid provided x and 
x' lie within a normal neighborhood^ [10' . However, this does not necessarily guarantee that a series representation will 
be convergent everywhere within this normal neighborhood. In fact, we will show that the series is only convergent 
within a smaller region, the size of which is given by the circle of convergence of the series. However, this does 
not preclude the use of the quasi-local expansion to calculate the Green function outside the circle of convergence 
(but within the normal neighborhood). As we will show, Pade resummation techniques, which have been extremely 
successful in other areas [TTl [T^ . are also effective in extending the series beyond its circle of convergence. 



'Electronic address; marc.casals@dcu.ie 
tElectronic address: sam.dolanOucd. iej 
Electronic address: adrian.ottewillOucd.ie' 
^Electronic address: barry.wardellOucd.ie 

^ More precisely, the Hadamard parametrix requires that x and x' lie within a causal domain - a convex normal neighbourhood with 
causality condition attached. This effectively requires that x and x' be connected by a unique non-spacelike geodesic which stays within 
the causal domain. However, as we expect the term normal neighbourhood to be more familiar to the reader, we will use it throughout 
this paper, with implied assumptions of convexity and a causality condition. 
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In this paper, we will focus in particular on calculating the Green function for two spherically symmetric spacetimes: 
Schwarzschild and Nariai. The Nariai spacetime |13|, I14j arises naturally from efforts to consider a simplified version 
of Schwarzschild [TJ [Kl dS]. It retains some of the key features of Schwarzschild (such as the presence of an unstable 
photon orbit and a similar effective radial potential which diminishes exponentially on one side) , but frequently yields 
more straightforward calculations. This makes it an ideal testing ground for new methods which are later to be 
applied to the more complicated Schwarzschild case. In the present work, we will use the line element of the static 
region of the Nariai spacetime (with cosmological constant A = 1 and Ricci scalar i? = 4) in the form 

ds^ = -(1 - p^)dt^ + (1 - p^y^dp^ + dnl, dnj = do^ + sin^ ed(t>'^. (1.3) 

where p G (—1, +1), t G (—00, +cxd), g [0, vr], G [0, 27r). In this form, it yields a wave equation with potential which 
is seen to closely resemble that of the Schwarzschild metric. 



ds^ ^ - i 1- —j dt^ + il - —j i 



dr^ + r^dni, (1.4) 



In Sec. [n] we use an adaptation of the Hadamard-WKB method developed by Anderson and Hu (T7] to efficiently 
calculate the coordinate series expansion of V{x, x') to very high order for both Nariai and Schwarzschild spacetimes. 



In Sec. Ill we use convergence tests to determine the radius of convergence of our series and show that, as expected, it 
lies within the convex normal neighborhood. We also estimate the local truncation error arising from truncating the 
series at a specific order. Using the method of Pade approximants, we show in Sec. |IV| how the domain of validity of 
the coordinate series can be extended beyond its radius of convergence to give an accurate representation of V{x, x') 
to within a small distance of the edge of the normal neighborhood. 



II. HADAMARD-WKB CALCULATION OF THE GREEN FUNCTION 



For the present quasilocal calculation, we need to consider the retarded Green function only for the points x and 
x' lying within a normal neighborhood. This allows us to express the retarded Green function in the Hadamard 
parametrix [HU [T5] , 

Gret {x, x) = 6*- {x, x) {U (x, x) 5 {<J (x, x)) — V (x, x) 6 {—(J (x, x'))} , (2.1) 

where 0_ {x,x') is analogous to the Heaviside step-function (i.e. 1 when x' is in the causal past of x, otherwise), 
5{a{x,x')) is the standard Dirac delta function, U {x^x') and V {x,x') are symmetric bi-scalars having the benefit 
that they are regular for x' x, and a{x,x') is the Synge [HI [HI [20] world function (i.e., half the square of the 
geodesic distance). The term involving U{x,x') is only non-zero for null connected points whereas the quasi-local 
self-force calculation which motivates us requires the Green function within the light-cone only. We will therefore 
only concern ourselves here with the calculation of the function V{x, x'). 

The fact that x and x' are close together suggests that an expansion of V{x,x') in powers of the separation of the 
points, 

oo 

V{x,x')^ ^ v,,^{r){t-t'f {cos ^-\y {r-r')\ (2.2) 

i,j,/c— 

where 7 is the angular separation of the points, may give a good representation of the function within the quasi- 
local region. Note that, as a result of the spherical symmetry of the spacetimes we will be considering, the expansion 
coefficients, Vijk{r), are only a function of the radial coordinate, r. Anderson and Hu [17 have developed a Hadamard- 
WKB method for calculating these coefficients. They applied their method to the Schwarzschild case and subsequently 
found the coefficients to 14*'' order using the Mathematica computer algebra system [2T]. In the present work, we 



adapt their method to allow for spacetimes of the Nariai form, (1.3). In particular, we consider a class of spacetimes 
of the general form 

ds^ = -f{r)dt^ + r\r)dr^ + g{r) [dO^ + sin^ Odcj)'^) , (2.3) 
where /(r) and g(r) are arbitrary functions of the radial coordinate, r, and previously Anderson and Hu had set 



g{r) = r , but allowed grr and gu to be independent functions of r. The form of the Nariai metric given in Eq. (1.3 1 
falls into the class (|2.3[), with /(r) = 1 — and g{r) = 1. For f{r) = 1 — ^ and g{r) — r^, this is the Schwarzschild 



metric of Eq. (1.4) 
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The method presented in this section differs from that of Ref. [T7] in the details of the WKB approach used, but 
otherwise remains very similar. Our alternative WKB approach, based on that of Refs. [22\ 123] proves extremely 
efficient when implemented in a computer algebra package. 

Following the prescription of Ref. [17 , the Hadamard parametrix for the real part of the Euclidean Green function 
(corresponding to the Euclidean metric arising from the change of coordinate t = it) is ^ 



-^r', x')] = ^ (^^P^ + V{x, x') H\a{x, x') |) + W{x, x')) , 
zn \ a{x, x') J 



(2.4) 



where U{x,x'), V{x,x') and W{x,x') are real-valued symmetric bi-scalars. 

Additionally, for the points x and x' separated farther apart in the time direction than in other directions. 



a{x,x') = ^-f{r){t^tr + 0[ix-xr] 



(2.5) 



so the logarithmic part of Eq. (2.4 1 is given by 



-V{x,x')ln{T-T'). 



(2.6) 



Therefore, in order to find V{x, x'), it is sufficient to find the coefficient of the logarithmic part of the Euclidean Green 
function. We do so by considering the fact that the Euclidean Green function also has the exact expression for the 
spacetimes of the form given in Eq. (2.3 1: 



G'e(— ir, x; —it 



1 /"OO ^ 

\x') = - dLJCOs[iu{T~T')]y^{2l+l)Pi{cos^)CM{r<)q^i{r>), (2.7) 



where p^ji and q^ji are solutions (normalised by C^ji) to the homogeneous radial equation for the scalar wave equation 
in the curved background (2.3) (and where r<, r> are the smaller/larger of r and r', respectively), along with the 
fact that 



du; cos Ht - r')] = ^^(r - r')^" log (r - r') 



(2n 
-1 
J2n) 



, (t-0'" log (T-r') + 



(2.8) 



where A is a low frequency cut-off justified by the fact that we will only need the log(r — r') term from the integral. 
We can therefore find V{x,x') as an expansion in powers of the time separation of the points by expressing the sum. 



^(2; + l)Pi(cos7)C^,p^/(r<)g^i(r>), 



(2.9) 



of Eq. (2.7) as an expansion in inverse powers of oj. This is achieved using a WKB-like method based on that 
of Refs. |22l [53]. Given the form (2.7) for the Euclidean Green function, the radial functions S{r) = Puji{r) and 
S{r) = quti{r) must both satisfy the homogeneous wave equation. 



^d's ^ 1 d ds 

dr"^ g dr 



dr 



T 



9 



'field 



s = o 



(2.10) 



where rrificid is the scalar field mass and ^ is the coupling to the scalar curvature, R. Next, given the Wronskian 
W{r) = C^i(p^iq'^i - q^ip'^i), its derivative is 



W = C^iiPujiqZi - qujip'U) 



1 

Tg 



{fg)'w 



(2.11) 



^ Note that this definition of the Green function differs from that of Ref. 1171 by a factor of 47r and the definition of V{x,x') differs by a 
further factor of 2. 
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and the Wronskian condition is therefore 

Cuji{Pu;iqii~ Qcuipii) = (2.12) 

19 

We now explicitly assume that r > r' and define the function 

B{r,r')^CM{r')qUr)- (2.13) 

Since the sum of Eq. (2.9 ) (and hence B(r, r')) is only needed as an expansion in powers of (r — r'), we expand B{r, r') 
about r' — r and (using Eq. (2.10) to replace second order and higher derivatives of i?(r,r') with expressions involving 
B{r,r') and dr'B{r,r')) find that 



B{r, r') = I3{r) + a{r){r' — r) + 



2(7? + X^) 



/3(r)-[ln(/5)]'a(r) 



(r' - rf 



(for r > r) (2.14) 



where 



and 



/?(r) = [i?(r,r')],,^,- = CMir)q^i{r), a(r) = [9,, i3(r, /)],,_,,- =C^(pL,(r)g^K'^) (2.15) 



,?(r)^--/g+(mLd + ei?)/5' (2.16) 
x2(r)^c.V + /.9f^+^') ■ (2.17) 



It will therefore suffice to calculate /?(r) and a(r). Furthermore, using Eq. (2.12) we can relate a{r) to the derivative 
of/3(r), 



a(r) 



V{r)9{ry 



(2.18) 



so it will, in f act, su ffice to find (3{r) and its derivative, /3'{r). 

Using Eqs. (2.10) and (2.12), it is immediate to see that P{r) must satisfy the nonlinear differential equati 



ion 



fg 



dr 



(2.19) 



The short distance behaviour of the Green function is determined by the high-cj and/or high-Z behaviour of the 
integrand of Eq. (2.7), so we seek to ex press (3{r) as an expansion in inverse powers of x- To keep track of this 
expansion we may replace x in Eq- (2.19 1 by x/e where e is a formal expansion parameter which we eventually set to 
1. Then, to balance at leading order we require 



4^3/2 



/3 



2x' 



We now write 



(2.20) 



(2.21) 



where /3o(^) = l/(2x(^))i insert this form for /3(r) in Eq. ( 2.19 1, and solve formally order by order in e to find a 
recursion relation for the (3n{r). On doing so and using Eq. (2.17) to eliminate (/+ i)^ in favour of and x^ we find 
that we can write 



PnV ) - 2^ 2n+2m+l 



m=0 



SO, for example, 



X 

A2,o(r-) ^ A2.i(r)a;2 ^2,2(0^^ ^ ^2,3('^)w^ ^ MAr) 
P2[r) = ^ + 7 + s + n + 13" 



(2.22) 

(2.23) 
(2.24) 
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The recursion relations for /3„(r) may then be re-expressed to allow us to recursively solve for the An^m{r)- Such a 
recursive calculation is ideally suited to implementation in a computer algebra system (CAS). Even on a computer, 
this recursive calculation becomes very long as n becomes large and, in fact, dominates the time required to calculate 
the series expansion of V{x, x') as a whole. For this reason, we have made available an example implementation 
in Mathematica, including precalculated results for several spacetimes of interest ^24. . This code calculates analytic 
results for A„ „i(r) on a moderate Linux workstation up to order n ^ 25, corresponding to a separation \x — x'\^'^ , in 
Schwarzschild and Nariai space-times in the order of 1 hour of CPU time. 

We also note at this point that knowledge of the A„ „i(r) (and their r derivative, which is straightforward to 
calculate) are all that is required to find the series expansion of /3'(r) and hence a(r). This can be seen by differentiating 
Eqs. (2.211 and (2.221 with respect to r to get 



with 

2n 



m=0 
2n+l 



,2n+2m+l 



{n + m+ -) 



1, 2xx'^«,m(r) 



(2.25) 



2n+2m+3 



- (n + TO - ]^)An,ra-l{r) 
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UJ 



2m 



X 



2n+2m- 



TT (2-26) 



where we have used Eq. (2.17) to write 2xx' = X^ifdY/ifd) + ^^((.9^)' ~ .'7^(/3)'/(/.9))i we use the convention 

^n.-l = ^n,2m+l — 0- 

With the An^mir) and their first derivatives calculated, we are faced with the sum over I in Eq. (2.9 1, where 
Eqs. ( |2.22| and '( |2.26| yield sums of the form 



E2(/+|)P,(cos7) 



^2n+2m+l ' 



1=0 



with 



An.ni (^) 



A'n.rnir) ~ {u + m + i)A„,™(r) ^^-{n + m- |)A„,„_i(r) (g^) 



^2y _ ifg)' „2 



fa 



(2.27) 



(2.28) 



Since we are considering the points x and x' to be close together, we can treat 7 as a small quantity and expand the 
Legendre polynomial in a Taylor series about 7 = 0, or, more conveniently, in powers of (cos 7— 1) about (cos 7— 1) = 0. 
It is straightforward to express each term in this series as a polynomial in even powers of {I + i): 



Pi(cos7) = 2F1 (-?, ; + 1; 1; (1 - cos7)/2) = ^ 

p=0 



((;+l)2„(l„l)2)...((^^l)2_(^„l)2) 

2P(p!)2 



(cos7-l)P (2.29) 



The calculation of the sum in Eq. (2. 7 1 therefore reduces to the calculation of sums of the form 



2D„,„(r)E 



^2n+2m+l ' 



(2.30) 



For fixed w and large I the summand behaves as " "'^ , and so only converges if p < n + ni. If p> n + m, we first 
split the summand as 



(/-hi)2p+iu;2™ 

■^2m+2n+l 



l)2(p-m-»)^2m / ^2g/f 



1-1/2 



(Z+ i)2(p-"-")c 

' (Jg)m+n+l/2 



,2m p-m-n 



E 

k=0 



(-l)'=(m + n+l/2)fe / Lu^g/f 



k\ 



-m—7i— 1/2 



P {-iy^{m + n+l/2)k f Lo-'g/f 

k=0 



kl 



(2.31) 



(2.32) 
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where {a)k = r(Q; + fc) /T{a) is the Pochhammer symbol and the sum correspond to the first {p — n — m) terms in the 
expansion of (1 + a;)-"-'"-i/2 about x = {uj^g/f)/{l + i)^ = 0. The terms outside the square brackets correspond 
to positive powers of oj and so contribute to the light cone singularity, not the tail term V{x, x') with which we are 
concerned in this paper. By contrast, the term in square brackets behaves as {I + |)~^ and so converges as Z — > cxj 
and will contribute to the tail term. We denote this term, with its prefactor as 



(Z+i)2p+ia;2m 



X 



2m+2n+l 



,2(p-n) 



Jp+l/2g2(ri+m)-p+l/2 ' 



m+n—p 



(l + xj 



-m—n — 1/2 



k=0 ^' / 



(2.33) 



and adopt the understanding that the sum vanishes if p < m + n. 
To proceed further, we use the Sommerfeld- Watson formula 



VF(;+i) = 5R - / dzF{z)itm{'Kz) ^ F 
1=0 I-W7 J -^0 



1 



(2.34) 



which is valid provided we can rotate the contour of integration for F(z) tan(7rz) from just above the real axis to the 
positive imaginary axis. Defining z = {f /gY^^ X/uj — l/^/a", where X = {I + 1/2), the sum (2.301 can then be written 
as the contour integral. 



^2(p-n) + l 
Jp+lg2(m+7i)-p 



dz 



(1 + 22)m+n+l/2 



2dz 



^2p+l 



_^ ^2nzuy/^ (1 - 22)™+»+l/2 



(2.35) 



Note that there is no need to include the rcgularization terms in the second integral as their contribution is manifestly 
imaginary and so will not contribute to the final answer. 



For p < m + n the first integral in Eq. (2.35 1 may be performed immediately as 

^2p+l pi 



(1 + z2)m+n+l/2 2(m + n - p - l/2)p+i ' 



(2.36) 



For p > m + n, we use the regularised integrand arising from Eq. (2.331 and temporarily introduce an ultraviolet 
cutoff l/e^ to get 



dx a;"+"-P-3/2 



(l + x) 



^ kl 

k=0 



-X 



(2.37) 



Integrating by parts p — m — n + 1 times, the regularisation subtraction terms ensure that boundary term go to zero 
in the limit e ^ and we are left with 



^_-^^p-rn-n+l (m + n + l/2)p_„,_„+i 



(1/2) 



p—7n—n-\-l 



dxa;-i/2(i + ^)-P-3/2 



(-l)P-'"-"+i7rp! 



r(m + n + l /2)T{p - 771 - n + 3/2) 



(2.38) 



where the last equality reflects that after these integrations by parts the remaining integral is finite in the limit e — > 0. 

The second integral in Eq. (2.351 is understood as a contour integral as illustrated in Fig. [l] The integrand is 
understood to be defined on the complex plane cut from z = — 1 to 1 and additionally possesses singularities at 
z — —1 and z = 1. To handle these in a fashion consistent with the Watson-Sommerfeld prescription we consider the 
contour in 3 parts: Ci running just below the cut from to 1 — e, C2 a semicircle of radius e about z — 1, and C3 
running along the real axis from z = 1 + e to 00. 

The integrand along C3 is manifestly imaginary and so gives zero contribution. Writing the integrand as 



G(z) 



where 



it is straightforward to see that 



3? 



L (1 



G{z) 



G{z) 



(1 - 2)™+"+l/2 



y2p+l 



(1 + g27r2(^V9//)(l + 2)™+"+l/2 



(1 — z)™+"+l/2 



dz 



E 

fc=0 



(-1)'= efc-™-"+l/2 



fc! 



I 1 ' 

m — n + i 



(2.39) 



(2.40) 



(2.41) 



7 



-> z 



c, 



FIG. 1: The second integral in Eq. (2.351 has a pole at z = I, so we split it into three parts: (1) An integral from to 1 
(2) An arc of radius e about z = 1 (3) An integral from 2 = 1 + e to oo. 



while integrating by parts m + n times 



(1 - z)™+"+l/2 



m+n-1 I ^ 

E ^ 

/£=0 



fc! 



k ^k—m—n-}~l/2 



dz 



(1 — 2)™+"+l/2 



G(z) 



*;=0 



fc! 



G('=)(1)(1-z)* 



(2.42) 



Adding the contributions from these components, it is clear that the e 

1 



m+n — 1 

E 

fe=0 



(-l)fc G(fe)(l) 
fc! fc — TO — n+i 



dz 



(1 — 2)ni+n+l/2 



G(z) 



divergences cancel and we are left with 
(-1) 



k=0 



fc! 



-G('=)(1)(1-z)* 



(2.43) 



where the subtraction terms in the integrand ensure the integral here is well-defined. 

While we can take this analysis further |26j . for our current purpose we note that we only need the expansion of 



Eq. (2.43 1 in terms of an inverse powers of uj as lu —>■ 00. From Eq. (2.40 1, it is immediate that the terms in the sum in 



(2.431 are exponentially small and so may be ignored for our purposes here. Indeed we may simultaneously increase 



the upper limit on the two sums in Eq. (2.431 without changing the result. Increasing it by 1 (or more), the integand 
increases from 0, peaks and then decreases to at 1 with the peak approaching as w — > 00. Standard techniques 
from statistical mechanics then dictate that the a; ^ 00 asymptotic form of the integral follows from the expanding 
the integrand, aside from the 'Planck factor', about z = and extending the upper limit to 00. Again, in doing so 
the contribution from the summation within the integrand give exponentially small contribution so that the powers 
of uj are determined simply by 



/ Series 

J (1 _|_ e^'^^'^v'?//) 



(1 — 22)m+ri+l/2 



(2.44) 



The coefficients in the series are known analytically, so to expand in inverse powers of lo we only need to compute 
integrals of the form 

■ dz (2.45) 



1 g27r2w^g// 



which have the exact solutions 



1 



r \b2n\ 



(2.46) 



where Bj^ is the A^-th Bernoulli number. This expression allows us to calculate the integrals very quickly. 

Applying this method for summation over I, Eq. (2.7 1 takes the form required by Eq. (|2.8|) so we now have 



the logarithmic part of the Euclidean Green function and therefore V{x, x') as the required power series in (t — 
t'), (cos 7 — 1) and (r — r'). Due to the length of the expressions involved, we have made available online [25 a 
Mathematica code implementing this algorithm, along with precalculated results for several spacetimes of interest 
including Schwarzschild, Nariai and Reissner-Nordstrom. 



8 



III. CONVERGENCE OF THE SERIES 

We have expressed V{x, x') as a power series in the separation of the points. This series will, in general, not be 
convergent for all point separations - the maximum point separation for which the series remains convergent will be 
given by its radius of convergence. In this section, we explore the radius of convergence of the series in the Nariai and 
Schwarzschild spacetimes and use this as an estimate on the region of validity of our series. 

For simplicity, we will consider points separated only in the time direction so we will have a power series in (t — t'), 



(3.1) 



ri=0 



where w„(r) is a real function of the radial coordinate, r only. We will also consider cases where the points are 
separated by a fixed amount in the spatial directions, or where the separation in other directions can be re-expressed 
in terms of a time separation, resulting in a similar power series in (t — t'), but with the coefficients, Vn{r), being 
different. This will give us sufficient insight without requiring overly complicated convergence tests. 



In the next section, we review some tests that will prove useful. In Sees. IIIB and III C we present the results of 
applying those tests in Nariai and Schwarzschild spacetimes, respectively. 



A. Tests for Estimating the Radius of Convergence 

1. Convergence Tests 



For the power series (3.1 1, there are two convergence tests which will be useful for estimating the radius of conver- 



gence. The first of these, the ratio test, gives an estimate of the radius of convergence, Atjic, 



At 



RC 



lim - 

n — *C30 



Vn+l 



(3.2) 



Although strictly speaking, the large n limit must be taken, in practice we have calculated enough terms in the series 
to get a good estimate of the limit by simply looking at the last two terms. The ratio test falls into difficulties, 
however, when one of the terms in the series is zero. Unfortunately, this occurs frequently for many cases of interest. 
It is possible to avoid this issue somewhat by considering non-adjacent terms in the series, i.e. by comparing terms 
of order n and n + to. 



(3.3) 



Using the ratio test in this way gives better estimates of the radius of convergence, although the results are still 
somewhat lacking. 

To get around this difficulty, we also use a second test, the root test, 



^tjic = lim sup 



1 



(3.4) 



which is well suited to power series. This gives us another estimate of the radius of convergence. Again, in practice 
the last calculated term in the series will give a good estimate of the limit. 

It may appear that only the (better behaved) root test is necessary for estimating the radius of convergence of the 
series. However, extra insight can be gained from including both tests. This is because for the power series (3.1 1, the 



ratio test typically gives values increasing in n while the root test gives values decreasing in n, effectively giving lower 
and upper bounds on the radius of convergence. 



2. Normal Neighborhood 



The Hadamard parametrix for the retarded Green function, (2.1 1, is only guaranteed to be valid provided x and x' 
are within a normal neighborhood (see footnote [T|. This arises from the fact that the Hadamard parametrix involves 
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Synge's world function, a{x,x'), which is only defined for the points x and x' separated by a unique geodesic. It is 
therefore plausible that the radius of convergence of our series could exactly correspond to the normal neighborhood 
size, ^NN- This turns out to not be the case, although it is still helpful to give consideration to t-^-^ as it should place 
an upper bound on the radius of convergence of the series. 

The normal neighborhood size will be given by the minimum time separation of the spacctime points such that they 
are connected by two geodesies. For typical cases of interest for self-force calculations there will be a particle following 
a time-like geodesic, so will be given by the minimum time taken by a null geodesic intersecting the particle's 
world-line twice. In typical black hole spacetimes, this geodesic will orbit the black hole once before re-intersecting 
the particle's world-line. 

Another case of interest is that of the points x and x' at constant spatial points, separated by a constant angle 
Acf). Initially, (i.e. when t = t'), the points will be separated only by spacelike geodesies. After sufficient time has 
passed for a null geodesic to travel between the points (going through and angle A(f>, they will be connected first by 
a null geodesic and subsequently by a sequence of unique time like geodesies. Since the geodesies are unique, ^nn will 
not be given by this first null geodesic time. Rather, tNN will be given by the time taken by the second null geodesic 
(passing through an angle 27r — A(j)). This subtle, but important, distinction will be clearly evident when we study 
specific cases in the next sections. 



3. Relative Truncation Error 



Knowledge of the radius of convergence alone does not give information about the accuracy of the series represen- 
tation of V{x,x'). The series is necessarily truncated after a finite number of terms, introducing a truncation error. 
As was done previously in Refs. [6l [28] the local fractional truncation error can be estimated by the ratio between the 
highest order term in the expansion ( O (At"), say) and the sum of all the terms up to that order, 

/ql M 

^ — \^ fa FT' (3-5) 
2^1=0 /ql W 

Refs. [HI [28] considered only the first two terms in the series when producing these estimates. Since we now have a 
vastly larger number of terms available, it is worthwhile considering these again to determine the accuracy of the high 
order series. 



B. Nariai Spacetime 

1. Normal Neighborhood 

Allowing only the time separation of the points to change, there are two cases of interest for the Nariai spacetime: 

1. The static particle which has a normal neighborhood determined by the minimum coordinate time taken by a 
null geodesic circling the origin (p = 0) before returning. This is the time taken by a null geodesic, starting at 
p — pi and returning to p' ^ pi, while passing through an angle A(j} = 2tt. 

2. Points at fixed radius, pi, separated by an angle, say tt/2. In this case, there is a null geodesic which goes 
through an angle A(j) — tt/2 when traveling between the points. However, there will not yet be any other 
geodesic connecting them, so this will not give tNN- Instead, it is the the next null geodesic, which goes through 
A(f) — 3n/2, that gives the normal neighborhood boundary. 

In both cases, the coordinate time taken by the null geodesic to travel between the points is given by [1] 

/l - pisech^(A(/)) + tanh(A</))A/l - /9fsech^(A0)\ 
tATAT = 2tanh-^(pi) -l-ln ^ . (3.6) 

\1 + pisech^(A(/)) - tanh(A(/))Wl - pjsech^ {A(j)) J 



2. Results 



For the Nariai spacetime, the ratio test suff'ers from difficulties arising from zeros of the terms in the series. This 
is because the i*^ term of the series (of order {t — t')^*) has 2i roots in p. In other words, the higher order the terms 
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FIG. 2: Radius of convergence as a function of the number of terms in the series for the Nariai spacetime with curvature 
coupling ^ — 1/6. The limit as n ^ oo will give the actual radius of convergence, but it appears that just using the terms up 
to n = 30 is giving a good estimate of this limit. The radius of convergence is estimated by the root test (blue dots) and ratio 
test (brown dots) and is compared against the normal neighborhood size (purple dashed line) calculated from considerations 
on null geodesies (see Sec. Ill A 2 and Sec. IIIBll. Note that plots for the ratio test were omitted in cases where it did not 
give meaningful results. 



considered, the more likely one of the coefficients is to be near a zero, and not give a useful estimate of the radius 
of convergence. We have therefore compared non-adjacent terms to avoid this issue as much as possible, by choosing 
m = n/2 in Eq. (3.3 1. Fortunately, the root test is much less affected by such issues and can be used without any 
adjustments. 

In Fig.|2] we fix the radial position of the points at p — p' — 1/2 and p — p' — 1/99 and plot the results of applying 
the root test (blue dots) and ratio test (brown dots) for cases ([Tl Static particle (left) and ([2| Points at fixed spatial 
points separated by an angle 7 = 7r/2 (right) as a function of the maximum order, rimax, of the terms of the series 
considered, i^nmax (t — t') It seems in both cases that the plot is limiting towards a constant value for the radius 

of convergence. In the case ([T]), we see that this gives a radius of convergence that is considerably smaller than the 
normal neighborhood size (purple dashed line). For the case ( [2|, we a gain see that the root test is limiting to a value 
for the radius of convergence. However, as discussed in Sec. Ill A 2 it is not the first null geodesic (lower dashed 
purple line) , but the second null geodesic (upper dashed purple line) that determines the normal neighborhood and 
places an upper bound on the radius of convergence of the series. 

In Fig.[3]we use the root test (blue dots) and ratio test (brown line)^ to investigate how the the radius of convergence 
of the series varies as a function of the radial position of the points, pi. Again, we look at both cases ([Tl (left) and ^ 
(right). As a reference, we compare to the normal neighborhood size (purple dashed line). We find that, regardless 
of the radial position of the points, the radius of convergence of the series is well within the normal neighborhood, 
by an almost constant amount. As before, in the case ^ of the points separated by an angle, we find that it is the 
second, not the first null geodesic that gives the normal neighborhood size. 

With knowledge of radius of convergence of the series established, it is also important to estimate the accuracy of 
the series within that radius. To that end, we plot in Fig. |4]the relative truncation error, (3.5 1, as a function of the 
time separation of the points at a fixed radius, p = p' = 1/2. We find that the 60"* order series is extremely accurate 



^ Note that the 'blips' in the ratio test are an artifact of the zeroes of the series coefficients used and are not to be taken to have any 
physical meaning. In fact, the 'blips' occur at different times when considering the series at different orders, so they should be ignored 
altogether. The ratio test plots should therefore only be fully trusted away from the 'blips'. Near the blips, it is clear that one could 
interpolate an approximate value, however we have not done so here as the plot is to be taken only as an indication of the radius of 
convergence. 
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FIG. 3: Estimates of the domain of validity (i.e. the radius of convergence) of the series expansion ofV{x,x') as a function 
of radial position m Nariai spacetime. Root test on O [(t — t')""] series given as blue dots, ratio test given by brown line (see 
footnote |3| and normal neighborhood estimate from null geodesies given by dashed purple lines. The left plot is for the case 
([l]), the static particle, right plot for the case (|2|, points separated by an angle of 7r/2. In both cases, the series is clearly 
divergent before the normal neighborhood boundary. This boundary is sometimes given by the second, rather than first null 
geodesic as can be seen in the plot on the right. In particular, this is the case for the points separated by an angle 7 = 7r/2 
since they are initially separated by only spacelike geodesies (see Sec. Ill B 1 1. 




FIG. 4: Relative truncation error in Nariai spacetime arising from truncating the series expansion for V{x, x') at order |a; — a;']^" 
(i.e. riniax ~ 30) for cases ([T]) the static particle (left) and ([2| points separated by an angle tt/2 (right). In both cases, the radial 
points are fixed at p = p' = 1/2. The series is extremely accurate until we get close to the radius of convergence (see FigB. 



to within a short distance of the radius of convergence of the series. 



C. Schwarzschild Spacetime 



1. Normal Neighborhood 



For a fixed spatial point at radius ri in the Schwarzschild spacetime, we would like to find the null geodesic that 
intersects it twice in the shortest time. This geodesic will orbit the black hole once before returning to ri. Clearly, the 
coordinate time ^nn for this orbit can only depend on ri. The periapsis radius, Vp, will be reached half way through 
the orbit. For the radially inward half of this motion, the geodesic equations can be rearranged to give 

r'" dr 
d<P=- , (3.7) 



! fl_2MA ^(1„ 2M) 



^""'^ dt . - r (3.8) 



For a given point ri, we numerically solve the first of these to find the periapsis radius, rp, then solve the second to 
give the normal neighborhood size, ^nn- 
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FIG. 5: Radius of convergence as a function of the number of terms considered. Root test (blue dots), ratio test (brown line - 
see footnote|3|, first null geodesic (purple dashed line). Left: Static particle at ri — lOOM Right: Static particle at ri = lOAf. 
The radius of convergence is given by the limit as the number of terms Wmax ^ oo. It is apparent that the plots are asymptoting 
to a value near, but slightly lower than the normal neighborhood size. Note that curves for the ratio test were omitted in cases 
where it did not give meaningful results. 
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FIG. 6: Radius of convergence as a function of radial position for a static point m Schwarzs child. Root test (blue dots), ratio 
test (brown line - see footnote [3|, first null geodesic (purple dashed line). 



2. Results 



The ratio test proves more stable for Schwarzschild than it was for Nariai. The series coefficients still have a large 
number of roots in r, but most are within r = 6Af and therefore don't have an effect for the physically interesting 
radii, r > 6M . 

Fig. |5] shows that the radius of convergence Atnc given by the root test is a decreasing function of the order of 
the term used, while that given by the ratio test is increasing. This effectively gives an upper and lower bound on 
the radius of convergence of the series. There is some 'noise' in the ratio test plot at lower radii (where that test is 
failing to give meaningful results), but we simply ignore this and omit the ratio test in this case. For the root test, 
the terms up to order {t — t')^"^ was used, while for the ratio test, adjacent terms in the series up order {t — t')^^ were 
compared. 

In Fig. |6]we apply the root (blue dots) and ratio (brown line - see footnote [s]) tests for the case of a static particle 
at a range of radii in Schwarzschild spacetime. Using the root test as an upper bound and the ratio test as a lower 
bound, it is clear that the radius of convergence is near, but likely slightly lower than the normal neighborhood size 
(purple dashed line). In this case, for the root test, the term of order {t — t')^"^ was used, while for the ratio test, 
terms of order {t — t')^'^ and {t — t')^^ were compared. 

In Fig. [7] we repeat for the case of the points separated on a circular timelike geodesic. The results are very similar 
to the static particle case and show the same features. 



IV. EXTENDING THE DOMAIN OF SERIES USING FADE APPROXIMANTS 

In the previous section it was shown that the circle of convergence of the series expansion of V{x, x') is smaller than 
the size of the normal neighborhood. This is not totally unexpected. We would expect the normal neighborhood size 
to place an upper limit on the radius of convergence, but we can not necessarily expect the radius of convergence to 
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FIG. 7: Radius of convergence as a function of radial position for a points separated along a circular geodesic in Schwarzschild. 
Root test (blue dots), ratio test (brown line - see footnote [sj , first null geodesic (purple dashed line). Note that curves for the 
ratio test were omitted in cases where it did not give meaningful results. 



be exactly the normal neighborhood size. However, since the Hadamard parametrix for the Green function is valid 
everywhere within the normal neighborhood, it is reasonable to hope that it would be possible to find an alternative 
series representation for V{x,x') which is valid in the region outside the circle of convergence of the original series, 
while remaining within the normal neighborhood. 

The radius of convergence found in the previous section locates the distance (in the complex plane) to the closest 
singularity of V{x,x'). However, that singularity could lie anywhere on the (complex) circle of convergence and will 
not necessarily be on the real line. In fact, given that the Green function is clearly not singular at the radius of 
convergence on the real line, it is clear that the singularity of V{x,x') does not lie on the real line. 

There are several techniques which can be employed to extend a series beyond its radius of convergence. Provided 
the circle of convergence does not constitute a natural boundary of the function, the method of analytic continuation 
can be used to find another series representation for V{x,x') valid outside the circle of convergence of the original 
series |29[l30j . This could then be applied iteratively to find series representations covering the entire range of interest 
of V{x,x'). Although analytic continuation should be capable of extending the series expansion of V{x,x'), there is 
an alternative method, the method of Pade approximants which yields impressive results with little effort. 

The method of Pade approximants [3TJ is frequently used to extend the series representation of a function beyond 
the radius of convergence of the series. It has been employed in the context of General Relativity data analysis with 
considerable success [TTl [12]. It is based on the idea of expressing the original series as a rational function (i.e. a ratio 
of two polynomials V{x,x') = R{x,x')/S{x,x')) and is closely related to the continued fraction representation of a 
function [31 . This captures the functional form of the singularities of the function on the circle of convergence of the 
original series. The Pade approximant, P^j{t — t') is defined as 

P^Cf t'] - ^"=0 ^"^^ ~ (A^] 

where Bq = 1 and the other {M + N + 1) terms are found by comparing to the first {M + + 1) terms of the original 
power series. The choice of M and N is arbitrary provided M + N < rt,„ax where rimax is the highest order term that 
has been computed for the original series. There are, however, choices for M and N which give the best results. In 
particular, the diagonal, , and sub-diagonal, Pl^_^i, Pade approximants yield optimal results. 
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FIG. 8: Comparing Pade approximant and Taylor series for 9{—a{x,x'))V(x,x') to 'exact' Green function from quasinormal 
mode sum in Nariai spacetime with curvature coupling ^ = 1/8 and ^ = 1/6. The Pade approximated V{x,x') (red line) 
is in excellent agreement with the quasinormal mode Green function (black dots) up to the normal neighborhood boundary, 
tNN ~ 6.56993 (top) and tNN ~ 4.9956 (bottom). 



A. Nariai 



The Green function in Nariai spacetime is known to be given exactly by a quasinormal mode sum [T] [33] at 
sufficiently late times. We can therefore use the Green function calculated from a quasinormal mode sum to determine 
the effectiveness of the Pade resummation. Figure [s] compares the quasinormal mode calculated Green function^ with 
both the original Taylor series representation and the Pade resummed series for V{x, x') for a range of cases. We 
use the Pade approximant P^^, computed from the 60"' order Taylor series. In each case, the series representation 
(blue dashed line) diverges near its radius of convergence, long before the normal neighborhood boundary is reached. 
The Pade resummed series (red line), however, remains valid much further and closely matches the quasinormal mode 
Green function (black dots) up to the point where the normal neighborhood boundary is reached. 

As was shown in Refs. ^[^, the Green function in Nariai spacetime is singular whenever the points are separated 
by a null geodesic. Furthermore, in Ref. [1] we have derived the functional form of these singularities and shown that 
they follow a four-fold pattern: 5{a), I/'kct, —5{a), — I/ttct, depending on the number of caustics the null geodesic has 
passed through (this was also previously shown by Ori [3S]). Within the normal neighborhood (where the Hadamard 
parametrix, (2.11, is valid), the 5{a) singularities (i.e. at exactly the null geodesic times) will be given by the term 
involving U{x,x ). However, a times other than the exact null geodesic times, the Green function will be given fully 
by V{x,x'). For this reason, we expect V{x,x') to reflect the singularities of the Green function near the normal 
neighborhood boundary 

The Pade approximant attempts to model the singularity of the function V(x, x') (which occurs at the null geodesic 
time) by representing it as a rational function, i.e. a ratio of two power series. By its nature, this will only faithfully 
reproduce singularities of integer order. In the Nariai case, however, the asymptotic form of the the singularities is 
known exactly near the singularity times, [I]- In cases where the points are separated by an angle 7 € (0,7r) (i.e. 
away from a caustic), the singularities are expected to have a l/{t — t' — tc) behavior and it is reasonable to expect 
the Pade approximant to reproduce the singularity well. When the points are not separated in the angular direction 
(i.e. at a caustic), however, the singularities behave like l/{t — t' — tc)"^/^ and we cannot reasonably expect the Pade 
approximant to accurately reflect this singularity without including a large number of terms in the denominator. 



* There are some caveats with how the quasinormal Green function was used. The fundamental mode (n = 0) Green function was used 
and a singularity time offset applied as described in Ref. [I]. In the case where two singularities are present, two sets of fundamental 
mode Green functions were used, each shifted by an appropriate singularity time offset and matched at an intermediate point. 
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FIG. 9: Relative Error in Improved vs Regular Fade Approximant. The relative error in the improved Pade approximant 
(orange) remains small closer to the singularity than the regular Pade approximant (green). 
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FIG. 10: Comparing Pade to Taylor series for Schwarzs child in the case of a static particle at r=10M. The Pade approximants 
Pii (red line) and PI5 (brown line) are likely to represent V{x,x') more accurately near the normal neighborhood boundary 
(at t — t' ^ 46. 2471 Af) than the regular Taylor series (blue dashed line). 

Given knowledge of the functional form of the singularity, however, it is possible to improve the accuracy of the 
Pade approximant further. For a singularity of the form 1/S{t), we first multiply the Taylor series by S{t). The result 
should then have either no singularity, or have a singularity which can be reasonably represented by a power series. 
The Pade approximant of this new series is then calculated and the result is divided by S{t) to give an improved 
Pade approximant. This yields an approximant which includes the exact form of the singularity and more closely 
matches the exact Green function near the singularity. In Fig. [9] we illustrate the improvement with an example 
case. We consider a static point in Nariai spacetime and compute the error in the Pade approximant relative to the 
quasinormal mode Green function (with n < 8) . The regular Pade approximant is shown in green while the improved 
Pade approximant is show in orange. The relative error remains small closer to the singularity for the improved Pade 
approximant case than for the the standard Pade approximant. Note that the error for early times arises from the 
failure of the quasinormal mode sum to converge and does not reflect error in the series approximations. 

B. Schwarzschild 

For the Schwarzschild case, there is no quasinormal mode sum with which to compare the Pade approximated series^. 
However, given the success in the Nariai case, we remain optimistic that Pade approximation will be successful for 
Schwarzschild. In an effort to estimate the effectiveness of the Pade approximant, we compare in Fig. [TO]the series 
expression for V{x,x') with two different Pade resummations, P|g and P||. The Pade approximant extends the 
validity beyond the radius of convergence of the series, but is less successful at reaching the normal neighborhood 
boundary {t — t' = Inn ~ 46.2471M) than in the Nariai case. 

The failure of the Pade approximant to reach the normal neighborhood boundary can be understood by the presence 
of extraneous singularities in the Pade approximant. The zeros of the denominator, S{t — t') = 0, give rise to 



^ A quasinormal mode sum could be computed for the Schwarzschild case, but would be augmented by a branch cut integral |T]. This 
calculation is in progress but has yet to be completed 
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singularities which occur at times earlier that the null geodesic time. It is possible that this problem could be 
reduced to a certain extent using the knowledge of the functional form of the singularities to compute an improved 
Pade approximant (as was successful in the Nariai case). However, to the authors knowledge, the structure of the 
singularities in Nariai spacetime is not yet known. While it may be possible to adapt the work of Ref. to find the 
asymptotic form of the singularities in Schwarzschild, without knowledge of the exact Green function we would not 
be able to dermine whether an improved Pade approximant would truly give an improvement. We therefore leave 
such considerations for later work. 



C. Convergence of the Pade Sequence 



The use of Pade approximants has shown remarkable success in improving the accuracy and domain of the series 
representation of V{x,x'). However, this improvement has not been quantified. There is no general way to determine 
whether the Pade approximant is truly approximating the correct function, V{x, x') or the domain in which it is valid 
|32j . In this subsection, we nonetheless attempt to gain some insight into the validity of the Pade approximants. 

The first issue to consider is the presence of extraneous poles in the Pade approximants. In Sec. |IV A| the Pade 
approximant was unable to exactly represent the l/{t — t' — tc)'^/'^ singularity at tc ~ 6.12 and instead represented it 
by three (real- valued) simple poles (at t — t' k, 6.522,6.854 and 9.488). This leads to the Pade approximant being a 
poor representation of the function near the poles. As was shown in Fig. [9] having exact knowledge of the singularity 
allows the calculation of an improved Pade approximant without extraneous singularities^. 

With extraneous poles dealt with, we consider the convergence of the Pade sequence of diagonal and sub-diagonal 
Pade approximants, 

P = {Pi P°, Pi Pi Pi Pi Pi, ■■■}, (4.2) 

with Pjv being the A^-th element of the sequence. The convergence of the Pade approximant sequence is determined 
by the behavior of the denominators, Sn for large N Provided S{x,x')n is not small, the Pade sequence will 
converge quickly toward the actual value of V{x,x'). When the first root of the denominator is at the null geodesic 
time (i.e. the normal neighborhood boundary), we can, therefore, be optimistic that the Pade sequence will remain 
convergent until this root is reached and the Pade approximants will accurately represent the function. 

To highlight the improvements made by using Pade approximants over regular Taylor series, we introduce the 
Taylor sequence (i.e. the sequence of partial sums of the series), T = {Tq, Ti,T2, ■ ■ ■ }, with A^-th element, 

N/2 

Tjv = ^«n(i-i')'"- (4.3) 

The two sequences P^r and Tn require approximately the same number of terms in the original Taylor series, so a 
direct comparison of their convergence will illustrate the improved convergence of the Pade approximants. 

In Fig. [TT] we plot the Pade sequence (blue line) and Taylor sequence (purple line) for the case of static points 
at p = 1/2 in the Nariai spacetime, with ^ = 1/8. For early times (eg. {t — t') = 2, it is clear that both Pade and 
Taylor sequences converge very quickly. At somewhat later times (eg. {t — t') — 3.3), both sequences appear to remain 
convergent, but the Pade sequence is clearly converging much faster than the Taylor sequence. Outside the radius of 
convergence of the Taylor series (eg. {t — t') — 5, 6.3), the Pade sequence is slower to converge, but appears to still do 
so. 



In Fig. 12 we again plot the Pade sequence, this time for the case of static points at r = lOAf in the Schwarzschild 
spacetime. As in the Nariai case, for early times (eg. {t — t') — lOM, it both Pade and Taylor sequences are converging 
very quickly. At slightly later times (eg. {t — t') = 20M, 27M) the convergence of the Pade sequence is better than 
the Taylor sequence. Outside the radius of convergence of the Taylor series (eg. {t — t') = 32M), the Pade sequence is 
slower to converge, but appears to still do so. The convergence of the series is slower in the Schwarzschild case than 
in the Nariai case. This is an indication that using more terms may yield a better result^. 



^ The improved Pade approximant has zeros in its denominator at t — f ' 2.64, 6.51, 6.59and8.73. The apparently extraneous singularity 
within the normal neighborhood (at t — t' x 2.64) does not cause any difficulty as the numerator also goes to zero at this point. 
In the Nariai case (with § = 1/8), the Taylor series for V{x, x') starts at order (t — t')^ and has been calculated to order {t — t')^'^, while 
for Schwarzschild it starts at (t — t')^ and has been calculated to order (t — t')^^. Since the Nariai series has several extra orders, it is 
reasonable to expect it to be a better approximation than the Schwarzschild series. 




FIG. 11: Convergence of the Taylor and Fade Sequences for the case of static points at p = 1/2 in the Nariai spacetime, 
with ^ = 1/8. Within the radius of convergence of the Taylor series, both Fade (blue line) and Taylor (purple line) sequences 
converge to the exact Green function (black dotted line) as calculated from a quasinormal mode sum with n < 6. The Fade 
sequence converges faster, particularly at larger times. Outside the radius of convergence of the Taylor series, only the Fade 
sequence is convergent. The top left plot is at a time {t — t') = 2, top right plot is at {t — t') = 3.3, bottom left at {t — t') — 5, 
bottom right at {t — t') — 6.3. 




FIG. 12: Convergence of the Taylor and Fade Sequences for the case of static points at r = lOM in the Schwarzschild spacetime. 
The top left plot is at a time {t - t') = lOAf, top right plot is at {t - t') = 20M, bottom left plot is at {t ~ t') = 27M, bottom 
right plot is at {t - t') = 32M. 



V. CONCLUSIONS 



In this paper we have presented an extension of the Hadamard-WKB method of Anderson and Hu [T^ to the Nariai 
spacetime. We have also demonstrated the use of an alternative WKB method [221 ES], which allows very high order 
terms in the Taylor series to be calculated efficiently (on a computer) for Schwarzschild, Nariai and other spherically 
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symmetric spacetimes. This allowed the series expansion of V{x,x') (appearing in the Hadamard parametrix of the 
Green function) to be computed to significantly higher order than was done previously in Refs. |171 121j. These high 
order expansions facilitated an investigation of the convergence properties of the series. We also demonstrated the 
huge benefit of Pade approximants to improving the domain and convergence of the series. 
This paper serves a dual purpose 

1. To discuss the calculation of the quasilocal Green function in Nariai spacetime, as required by Ref. 

2. To investigate the potential for applying the same techniques in the Schwarzschild spacetime, with the goal of 
computing an accurate quasilocal Green function for use in a matched expansion calculation of the self- force. 

We have found that, using Fade approximants, it is possible compute the quasilocal Green function in Nariai 
spacetime to high accuracy to within a short distance of the normal neighborhood boundary. Even without the use of 
Fade approximants, the Taylor approximated series gives good accuracy within a large part of the quasilocal region. 
This gives confidence in their use for matched expansion calculations in Ref. [1]. 

With regard to the Schwarzschild case, we find that the quasilocal calculation of the Green function is in good 
standing and should be usable in matched expansion calculations once techniques for computing the 'distant-past' 
Green function have been fully developed. Both Fade and Taylor sequences remain convergent within a large part 
of the normal neighborhood. The use of Fade approximants is not quite as successful as for the Nariai case, but 
we remain optimistic that knowledge of the structure of the singularities in Schwarzschild may allow for the use of 
improved Fade approximants as discussed in Sec. |IVB[ 

The orders of the series calculated for Nariai (60-th) and Schwarzschild (52-nd) were the maximum possible within 
a reasonable time (~ 1 day on a modern Linux desktop). Although these are considerably high order series, one may 
still wonder whether they are sufficiently high for matched expansion calculations. It is clear from Sec. |IV C| that the 
higher order terms only have a significant contribution near the radius of convergence (for the Taylor series) or normal 
neighborhood boundary (for the Fade approximant) . As the Hadamard parametrix is only valid within the normal 
neighborhood, we consider the fact that the Fade approximant is accurate to within a short distance of the normal 
neighborhood boundary to be confirmation that the series has been calculated to sufficiently high order (in particular 
for the Nariai case). Additionally, within the matching region used in Ref. [T] for Nariai, we clearly have computed a 
sufficient number of coefficients to give the Green function to high accuracy. We can be optimistic that this is also 
the case for Schwarzschild: the quasilocal series is accurate long after the time when the quasinormal mode sum is 
expected to be convergent {t — t' — 2r* « 12.77Af, for the case of a static particle at r = lOM considered here). 

Our analysis has remained focused primarily on the case of one dimensional series. This was done for reasons of 
simplicity and clarity. For multi-dimensional series, one could re-express each of the coordinates in terms of a single 
parameter, as was done in Sec. |III| for the case of a circular geodesic in Schwarzschild. Alternatively, one could make 
use of the extension of the Fade approximant to double and higher dimensional power series as developed by Chisholm 

[ISIST]. 

VI. ACKNOWLEDGMENTS 

MC is grateful to the Department of Fhysics and Astronomy of the University of Mississippi for its hospitality 
during the preparation of this paper. MC was partially funded by Fundagao para a Ciencia e Tecnologia (FCT) - 
Portugal through project FTDC/FIS/64175/2006. MC, BW and SD are supported by the Irish Research Council for 
Science, Engineering and Technology, funded by the National Development Flan. 



[1] M. Casals, S. Dolan, A. C. Ottewill, and B. Wardell (2009), arXiv:0903.0395. 
[2] T. C. Quinn and R. M. Wald, Phys. Rev. D56, 3381 (1997), gr-qc/9610053. 
[3] T. C. Quinn, Phys. Rev. D62, 064029 (2000), gr-qc/0005030. 

[4] Y. Mino, M. Sasaki, and T. Tanaka, Phys. Rev. D55, 3457 (1997), gr-qc/9606018. 
[5] B. S. DeWitt and R. W. Brehme, Ann. Phys. 9, 220 (1960). 
[6] A. C. OttewiU and B. Wardell, Phys. Rev. D7T, 104002 (2008), 0711.2469. 
[7] A. C. OttewiU and B. Wardell, Phys. Rev. DT9, 024013 (2009), 0810.1961. 
[8] E. Poisson, Living Rev. Relativity 7, 6 (2004), gr-qc/0306052. 
[9] S. Detweiler, Class. Quantum Grav. 22, S681 (2005), gr-qc/0501004. 
[10] F. G. Friedlander, The Wave Equation on a Curved Space-time (Cambridge University Press, Cambridge, 1975), ISBN 
978-0521205672. 



19 



T. Damour, B. R. Iyer, and B. S. Sathyaprakash, Phys. Rev. D 57, 885 (1998). 

E. K. Porter and B. S. Sathyaprakash, Phys. Rev. D 71, 24017 (2005), gr-qc/0406038. 

H. Nariai, Sci. Rep. Tohoku Univ. 34, 160 (1950). 

H. Nariai, Sci. Rep. Tohoku Univ. 35, 62 (1951). 

V. Cardoso and J. P. S. Lemos, Phys. Rev. D 67, 084020 (2003), gr-qc/0301078. 
S. Zerbini and L. Vanzo, Phys. Rev. D 70, 044030 (2004), hep-th/0402103. 
P. R. Anderson and B. L. Hu, Phys. Rev. D69, 064039 (2004), gr-qc/0308034. 

J. Hadamard, Lectures on Cauchy's Problem in Linear Partial Differential Equations (Dover Publications, 1923), ISBN 
978-0486495491. 

J. L. Synge, Relativity: The General Theory (North-HoUand, Amsterdam, 1960), ISBN 978-0720400663. 
B. S. DeWitt, Dynamical theory of groups and fields (Gordon and Breach, New York, 1965). 
P. R. Anderson, A. Eftekharzadeh, and B. L. Hu, Phys. Rev. D73, 064023 (2006), gr-qc/0507067. 
K. W. Howard, Phys. Rev. D30, 2532 (1984). 

E. Winstanley and P. M. Y oung, Phys. Rev. D77, 024 008 (2008 ), 0708.3820. 



http: //www. barrywardell .net /research/ code/hadamard-wkb 

G. Watson, Proc. Roy. Soc. Lond. A 95, 83 (1918) 

A. C. OttewiU, E. Winstanley, and P. M. Young (2009), in preparation. 

I. Gradshteyn and I. Ryzhik, Table of Integrals, Series, and Products (Academic Press, 2007). 

W. G. Anderson and A. G. Wiseman, Class. Quantum Grav. 22, S783 (2005), gr-qc/0506136. 

P. Morse and H. Feshbach, Methods of Theoretical Physics (McGraw-Hill Book Company, 1953). 

E. T. Whittaker and G. Watson, A Course of Modern Analysis (Cambridge University Press, 1927), 4th ed. 

C. M. Bender and S. A. Orszag, Advanced Mathematical Methods for Scientists and Engineers (Springer, 1999). 

W. H. Press, S. A. Teukolsky, W. T. Betterling, and B. P. Flannery, Numerical Recipes (Cambridge University Press, 

2007), 3rd ed. 

H. R. Beyer, Communications in Mathematical Physics 204, 397 (1999). 

B. Kay, M. Radzikowski, and R. Wald, Communications in Mathematical Physics 183, 533 (1997). 

A. Ori (2008), The four-fold structure of the singular part of the Green's function beyond the caustics was found by Ori 
some time ago. Ori conducted measurements in analog acoustic system, which seem to verify the theoretical prediction (A. 
Ori, Private Communication, 2008). 

J. S. R. Chisholm, Mathematics of Computation 27, 841 (1973). 

J. S. R. Chisholm and J. McEwan, Proc. R. Soc. Lond. A 336, 421 (1974). 



